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ABSTRACT 



Aims. We study the dependence of galaxy clustering on luminosity and stellar mass at redshifts z ~ [0.2 - 1], using the first lOK redshifts from the 
zCOSMOS spectroscopic survey of the COSMOS field. 
' Methods. We measure the redshift-space correlation functions ^(rp,n) and ^(s) and the projected function, Wp(rp) for sub-samples covering 
different luminosity, mass and redshift ranges. We explore and quantify in detail the observational selection biases due to the flux-limited nature 
of the survey, using ensembles of realistic semi-analytic mock samples built from the Millennium simulation. We use the same mock data sets to 
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^ carefully check our covariance and error estimate techniques, comparing the performances of methods based on the scatter in the mocks and on 
I ^ I bootstrapping schemes. We finally compare our measurements to the cosmological model predictions from the mock surveys. 

Results. At odds with other measurements at similar redshift and in the local Universe, we find a weak dependence of galaxy clustering on 
I luminosity in all three redshift bins explored. A mild dependence on stellar mass is instead observed, in particular on small scales, which becomes 
^ particularly evident in the central redshift bin (0.5 < z < 0.8), where Wp(rp) shows strong excess power on scales > 1 h"' Mpc. This is reflected 
' in the shape of the full ^{ri„n) that we interpret as produced by large-scale structure dominating the survey volume and extending preferentially 
in direction perpendicular to the line-of-sight. Comparing to z ~ measurements, we do not see any significant evolution with redshift of the 
amplitude of clustering for bright and/or massive galaxies. 
I Conclusions. This is consistent with previous results and the standard picture in which the bias evolves more rapidly for the most massive halos, 
• which in turn host the highest-stellar-mass galaxies. At the same time, however, the clustering measured in the zCOSMOS lOK data at 0.5 < z < 1 
for galaxies with log(M/MQ) > 10 is only marginally consistent with the predictions from the mock surveys. On scales larger than ~ 2 h ' Mpc, 
the observed clustering amplitude is compatible only with ~ 1% of the mocks. Thus, if the power spectrum of matter is ACDM with standard 
, normalization and the bias has no "unnatural" scale-dependence, this result indicates that COSMOS has picked up a particularly rare, ~ 2 - 3cr 
' positive fluctuation in a volume of ~ 10* h"' Mpc^. These findings underline the need for larger surveys of the z ~ 1 Universe to appropriately 
^ ' characterize the level of structure at this epoch. 

Key words. Cosmology: observations - Cosmology: large-scale structure of Universe - Surveys - Galaxies: evolution 

1. Introduction hosting halo is expected to play a significant role in the defini- 
tion of the visible properties of the galaxy, as the total mass in 

In the canonical scenario of galaxy formation, galaxies are ^^^^^^ -^^ luminosity, color, star formation rate and, pos- 

thought to form through the cooling of b aryonic gas within ex- ^j^jy morphology 
tended dark matter halos (IWhite & ReesL lT978). The mass of the 

Since it is the baryons that form the visible fabric of the 

Send offprint requests to: B. Meneux bmeneux@mpe.mpg.de Universe, a major challenge in testing the galaxy formation 
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erties one is actually constraining the relationship between the 
dark and luminous components of galaxies. 

Measurements of first moments, as the luminosity function 
or the stellar mass function, provide a way to understand how 
these are related to the total halo mass fu nctions, which can 
be o btained from analytic predict ions (e.g. Press & SchechteA 
11974 ) or n-body simulations (e.g. IWarren et al.L l2006l) . Similar 
investigations can be made on the se cond moment, i.e. th e 
two-point correlation function (e.g. ISpringel et al.L l2006l) . 
Studies of galaxy clustering in large local surveys have shown 
how clustering at z ~ does depend significant ly on sev- 
eral s p ecific propert ie s. The se include lu minosity (Hamilton 



19881: llovino et al.L 119931: iMaurogordato & Lachieze-Revt 
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dWill meret al., 1998; Norb erg et al 
morphology (iDayis & Gellei 



2002) 



20001: iNorberg et al 
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[1976 



^ , Giovanel U et al „ 

s tellar mass tLi et al.> .2006F and 



ogy (lUay: 

19861 iGuzzoetaLl Il997h, - - -- - 
environment (lAbbas&She"fliLl2006h . 

In recent years it has become possible to extend these 
investigations to high redshift, obtaining first ind icative re 



suits o n how thes e dependences evolve with t i me jCoil et al 
2006[ JPpllo et al.L 120061: JPhleps et al.', '20061 iMeneux et al" 
2006t iDaddi et al] l2003t [Meneux et al., 200^. The VIMOS- 
VLT D eep Survey (VVD S) (IPolloet all 120061) and the DEEP2 
survey (ICoil et al.L l2006l) in particular, have provided new in- 
sights on the way galaxi es of differe nt luminosity cluster at 
z ~ 1 . More specifically, i PoUo et al] ([2006) have shown that 
at these epochs galaxies akeady show a luminosity segregation, 
with more luminous galaxies being more clustered than faint ob- 
jects. At the same time, however, a significant steepening with 
luminosity of the shape of their two-point correlation function 
for separations < 1 - 2 h"' Mpc, is observed. This behaviour 
is at variance with that at z ~ 0. A similar trend h as been ob- 
serve d at the same redshift by the DEEP2 sur vey dCoil et al.L 
|2006|) . Complementarity, Meneux et al.l (l2008l) have shown a 
positive trend of clustering with stellar mass also at z ~ 1 , with a 
clear evidence for a stronger evolution of the bias factor for the 
most massive galaxies (see also Brown et al., 2008; Wak e et all 
120081) . 

The interpretation of the evolution in shape and amplitude 
of Wp{rp) with respect to luminosity and redshift is particularly 
interesting in the context of the halo model for galaxy forma- 
tion. In this framework, the observed shape of ^(r) (or Wpirp)) is 
interpreted as composed by the sum of two components: (a) the 
1-halo term, which dominates on small scales (< 1 -2 h"' Mpc at 
the current epoch), where correlations are dominated by pairs of 
galaxies living within the same dark-matter halo (i.e. in a group 
or cluster); (b) the 2-halo term on large scales, which is char- 
acterized by pairs of galaxies occ upying different dark-matter 
halos (see ICoorav & ShethI (l2002 l) for a review). IZheng et"aLl 
(l2007h have modelled the luminosity-dependent Wp{rp) from 
both the DEEP2 (at z ~ 1) and SDSS (at z ~ 0) surveys, within 
such Halo Occupation Distribution (HOD) framework. In this 
way they establish evolutionary connections between galaxies 
and dark-matter halos at these two epochs, providing a self- 
consistent scenario in which the growth of the stellar mass de- 
pends on the halo mass. Similar results are obtained more re- 
cently in a combine d analysis of the VVDS-Deep and SDSS data 
(lAbbas et al.Ll2009h . 

In this paper we use the first 10,000 redshifts from the zCOS- 
MOS redshift survey (the "lOK sample") to further explore 
these high-redshift trends of clustering with luminosity and mass 
based on a new, independent sample. Although shallower than 



VVDS-Deep and DEEP2 {Ub < 22.5 vs. 24 and 23.5, respec- 
tively), zCOSMOS covers a significantly larger area and samples 
a volume of ~ 3 x 10^ h ' Mpc to redshift z=1.2. This should 
hopefully help reducin g the effect of cosmic variance (still stron g 
for samples this size ( GarilU et al.L 120081; IStringer et al.Ll2009h ). 
while providing a better sampling of the high-end tail of the lu- 
minosity and mass functions. However, one main result from this 
analysis will be the explicit demonstration of how strong cos- 
mic variance still is within volumes of the Universe this size. 
The clustering properties of the zCOSMOS sample in the vol- 
ume contained within the redshift range 0.4 - 1 seem to lie at 
the extreme high end of the distribution of fluctuations on these 
scales, as it was in fact alread y suggested by the angu lar clus- 
tering of the COSMOS data dMcCracken etalL |2007|) . As we 
shall see, these results an d those presen ted in th e zCOSMOS se- 
ries of clustering papers ([Porciani et"a l., 2009; de la Torre et al.L 
l2009l:lAbbas iraill2009h indicate how cautious one should be in 
drawing far-reaching conclusions from the modelling of current 
clustering results from deep galaxy surveys. 

A significant part of this paper is dedicated to discussing in 
detail these cosmic-variance effects, together with the impact 
of incompleteness on the derived results. This is particularly 
important when constructing mass-limited sub-samples from a 
magnitude-limited survey, which introduces a mass incomplete- 
ness that depends on redshift and stellar mass. The intrinsic scat- 
ter in the galaxy mass-luminosity relation determines a progres- 
sive loss of faint galaxies with high mass-to-light ratio. We study 
in detail the eff'ect of this incompleteness on the measured clus- 
tering both using the data themselves and mock samples built 
from the Millennium simulation. At the same time, we explore 
in quite some detail our ability to characterize measurement er- 
rors and the covariance matrix of our data, comparing estimates 
from the mock samples to those from bootstrap resamplings of 
the data themselves. 

The paper is organised as follows. In Sect. 2 and 3 we de- 
scribe the zCOSMOS survey and the simulated mock samples 
used in the analysis, while in Sect. 4 we describe the selection 
of luminosity- and mass-limited subsamples, discussing exten- 
sively the incompleteness related to this operation; in Sect. 5 we 
describe our clustering estimators, while in Sect. 6 we discuss 
in detail the observational biases and selection effects, how we 
account for them and what is their effect on the measured quan- 
tities; in Sect. 7 we explore and discuss the error budget and 
how to estimate the covariance properties of our measurements; 
in Sect. 8 and 9 we present our measurements of clustering as a 
function of luminosity and mass, respectively, while in Sect. 10 
we compare these results with those from other surveys and with 
simple model predictions; finally, in Sect. 11 we place these find- 
ings in a broader context and discuss future developments. 

Throughout the paper we adopt a cosmology with Q,„ = 
0.25, Qa = 0.75. When needed, we also adopt a value erg = 0.9 
for the normalization of the matter power spectrum; this is cho- 
sen for consistency with the Millennium simulation, also used 
for comparison to model predictions. The Hubble constant is pa- 
rameterized via h - //o/lOO to ease comparison with previous 
works. Stellar masses are quoted in unit of h - 1. All length 
values are quoted in comoving coordinates. 

2. The zCOSMOS survey data 

The zCOSMOS survey dLillv et al.L l2007h is performed with 
the VIMO S multi-object spectr ograph at the ESO Very Large 
Telescope (iLe Fevre et al.L[2()03h . 600 hours of observation have 
been allocated to this program. These are invested to mea- 
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Fig. 1. Distribution on the sky of the ~ 10,000 galaxies with 
measured redshift (crosses) forming the zCOSMOS "lOK" sam- 
ple. The large blue dots mark the centres of independent VIMOS 
pointings, each including four quadrants on the sky (as described 
by the red solid lines). 



sure sp ectra for galaxies in the COSMOS field dScoville et al.L 
l2007 a^). targeting: (a) -20000 galaxies brighter than / < 22.5 
(zCOSMOS Bright); -10000 sources at redshift 1.4 < z < 3.0 
pre-selected using color-color criteria (Lilly et al.L l2007l) (zCOS- 
MOS Faint). So far, the survey has observed about half of the to- 
tal "Bright" sample. This is the so called "lOK" sample used for 
the analysis presented in this paper, and is based on the observa- 
tions of 83 VIMOS poi ntings over 44 distinct telescope positions 
on the sky (iLillv et al.L 12009 ). These are shown in Fig.[Tl where 
the footprint of VIMOS (4 quadrants of ~ 7 x 8 arcmin- sepa- 
rated by a cross about 2 arcmin^ wide) is evident. About every 
3'''' galaxy has been observed in the field. The final "20K" zCOS- 
MOS sample will be twice larger, reaching a sampling around 
60-70%. The correction of the complex angular selection func- 
tion will be discussed later in the context of our galaxy clustering 
measurement. 

Observations are performed using the medium resolution 
RED grism, corresponding to R ~ 600 and covering the spectral 
range 5550 - 9650A. The average error on the redshift measure- 
ments has been estimated from the repea ted observations o f 632 
galaxies, and is found to be ~100 km s ' (iLillv et al.Ll2009l) . This 
roughly corresponds to a radial distance error of 1 h ' Mpc. The 
reduction of the data to the redshift assignment was carried out 
independently at two institutes before a reconciliation process to 
solve discrepancies. The quality of each measured redshift was 
then quantified via a quality flag that p rovides us with a confi- 
dence level (see ILillv etal] (l2007l |2009|) for definition). For the 
present work, we only use redshifts with flags 1.5 to 4.5 and 9.3 
to 9.5, corresponding to confidence levels greater than 98%. 

The zCOSMOS survey benefits of the la rge multi- 
wavelength coverage of the COSMOS field dCapak et al., 2007), 
that w ith the latest ad ditions now comprises 30 photometric 
bands jllbert et al.L 120091) extending well into the infrared. These 



include in particular accurate /T-band and Spitzer-IRAC pho- 
tometry over the whole area, which have allowed us to de- 
rive relevant physical properties as rest-fra me luminos i ty and 
stellar mass witli unprecedented accuracy dZucca et al.1 120091; 
iBolzonella et"!!]. 120091: IPozzetti et al.Ll2009l) . 

3. Mock survey catalogues 

In this paper we shall make intense use of mock surveys con - 
structed from the Millennium simulation (ISpringel et al.Ll2005h . 
This is done (a) to understand the effect of our selection criteria 
on the measured quantities (Sec. |6.3l i: (b) to estimate the mea- 
surement errors and covariance of the data (Sec.|7]i. 

We use two sets of l ight c ones, constructed a s expla ined 
in Kitzb ichler & White! (l2007l) and iBlaizot et al.l (12005) by 
combining dark-matter halo trees from the Millennium run 
to the Munich semi-anal ytical model of galaxy formation 
(iDe Lucia & B laizot*. '2007^. The two sets contain 24 1.4 x 1.4 
deg2 mocks built by Kitzbichler & White (2007) and 40 1 x 1 
deg^ mocks built by D e Lucia & BlaizotI (12007) . that we shall 
name KW24 and DLB40, respectively. The main difiference be- 
tween the two sets, in addition to the different survey area, is 
that the DLB40 set contains all galaxies irrespective to any cri- 
teria down to the simulation limit that corresponds roughly to 
10** Mq, up to redshift z - 1.7, whereas the KW24 set con- 
tains only galaxies brighter than / < 22.5. This implies that the 
DLB40 set allows us to select in stellar mass down to very low 
masses and test selection effects. The observing strategy of the 
zCOSMOS lOK sample was only applyied to the KW24 set, al- 
lowing us to do a carefull error analysis of our measurements. 

The Millennium run contains = 2160^ particles of mass 
8.6 X 10*^ h"' Mo in a cubic box of size 500 h~' Mpc. The 
simulation was built with a ACDM cosmological model with 

= 0.25, Qa = 0.75, erg = 0.9 and = 73 km s"' Mpc"'. 

4. Luminosity- and IVIass-selected subsamples 

4. 1 . Luminosity selection 

Absolute magnitu des have been derived for the IQK g alaxies us- 
ing the code ALF dllbert et all 120051: IZucca etall 12009). which 
is based on fitting a Spectral Energy Distribution to the observed 
multi-band photometry. There are various sources of uncertain- 
ties to take into account (errors on apparent magnitudes, num- 
ber of available photometric bands, method used, . . . ). A direct 
comparison with absol ute magnit udes derived with the indepen- 
dent code ZEBRA (Fel dmann et al.. 2006) shows consistent es- 
timates with a small dispersion of cr ~ 0.05 magnitudes, in par- 
ticular in the B band. This can resonably be considered as the 
typical error on our absolute magnitudes. 

For our analysis, the goal is to define luminosity-limited 
samples that are as close as possible to truly volume-limited 
samples, i.e. with a constant number density. This should be pos- 
sibly done within a few independent redshift ranges. The size of 
the redshift slices in which to split the sample has to be chosen 
as a compromise between two aspects: (a) we want it to be large 
enough as to have sufficient statistics and provide a good mea- 
surement of clustering; however, (b) we do not want it to be too 
large, as to avoid significant evolution within each redshift bin. 

However, we know that all through the overall redshift range 
covered by the zCOSMOS survey (0.2 < z < 1.1) the luminos- 
ity of galaxies evolves, with a clear ch ange in the c haracteristic 
parameters of the luminosity function dllbert et al.L (20051) . This 
evolution does depend on the morphological/spectral type of the 
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Fig. 2. Selection boundaries of the different sub-samples of the zCOSMOS lOK survey used in this paper Left: luminosity-redshift 
selection, which accounts for the average luminosity evolution of galaxies; Right: mass-redshift selection. 



galaxy considered. To be able to select a nearly volume-limited 
sample within a given redshift interval, we need to take the cor- 
responding evolution into account. This can be realistically done 
only in a statistical way, looking at the population-averaged evo- 
lution of the global luminosity function. 

We have therefore considered t he observed lumin osity func- 
tion measured from the same data jZucca et al.Ll2009l) and mod- 
elled its change with redshift as a pure luminosity evolution (i.e. 
keeping a constant slope a and normalisation factor <t)*), which 
is a fair description of the observed behaviour We find that the 
characteristic absolute magnitude M*{z) evolves with redshift as 

M*(z)^Ml+Az , (1) 

where A 1 . In the companion paper, Ide la Torre et"ai] (l2009l) 

split the zCOSMOS galaxy samples in 3 morphological classes. 
They observe different luminosity evolutions for elliptical, spi- 
ral and irregular galaxies, with A varying from ~ -0.7 for to 

1.2 but wit h large uncerta i nties m aking A - -\ compatible 

for all classes. iPorciani et alj (l2009l) reach similar conclusions 
when dividing the zCOSMOS lOK sample into 3 color classes. 

We therefore define our luminosity-limited samples by an ef- 
fective absolute magnitude cut at z = 0, Mb,cui and including all 
galaxies with Mb{z) - 5log(h) < Mb^cm - Z- The resulting selec- 
tion loci for different values of Mb^cui are plotted over the data 
in the luminosity-redshift plane in the left panel of Fig. |2] As 
evident from the figure, the faintest allowed threshold Mb.cui de- 
pends on the redshift range considered, i.e. z=[0.2-0.5], z=[0.5- 
0.8] and z=[0. 8-1.0]. The details of the resulting samples are 
described in Table [1] 

4.2. Mass selection 

Stellar mass has become a quantity routinely measured in re- 
cent years, thanks to surveys with mu lti- wavelength pho tome- 
try, extending to the near-infrared (e.g. iRettura et all 12006) . al- 



Table 1. Properties of the luminosity-selected samples 



Sample 


Redshift 


Mean 


Mb,cu, 


Number of 




range 


redshift 


(z=0) 


galaxies 


Ll.l 


0.2-0.5 


0.37 


-18.00 


1892 


L1.2 


0.2-0.5 


0.37 


-18.50 


1311 


L1.3 


0.2-0.5 


0.37 


-19.00 


811 


L1.4 


0.2-0.5 


0.37 


-19.50 


469 


L2.1 


0.5-0.8 


0.67 


-19.00 


1848 


L2.2 


0.5-0.8 


0.67 


-19.50 


1025 


L2.3 


0.5-0.8 


0.67 


-20.00 


441 


L3.1 


0.8-1.0 


0.91 


-19.50 


971 


L3.2 


0.8-1.0 


0.91 


-20.00 


447 



Table 2. Properties of the mass-selected samples 



Sample 


Redshift 


Mean 


log(M/Me) 


Number of 




range 


redshift 


range 


median 


galaxies 


Ml.l 


0.2-0.5 


0.36 


> 9.0 


9.80 


2159 


Ml. 2 


0.2-0.5 


0.37 


> 9.5 


10.09 


1445 


Ml. 3 


0.2-0.5 


0.36 


> 10.0 


10.36 


827 


M1.4 


0.2-0.5 


0.37 


> 10.5 


10.66 


275 


M2.1 


0.5-0.8 


0.66 


> 9.0 


9.97 


2831 


M2.2 


0.5-0.8 


0.66 


> 9.5 


10.12 


2276 


M2.3 


0.5-0.8 


0.67 


> 10.0 


10.38 


1366 


M2.4 


0.5-0.8 


0.67 


> 10.5 


10.68 


477 


M3.1 


0.8-1.0 


0.90 


> 10.0 


10.46 


755 


M3.2 


0.8-1.0 


0.90 


> 10.5 


10.73 


344 



though some uncertaintie s related to the detai led modelling of 
stellar evolution remain (iPozzetti et all l2007l) . This has made 
studies of clustering as a function of stellar mass possible for 
large statistical samples. We used stellar masses estimated by 
fitting the Spectral Energy Distribution (SED), as sampled by 
the large multi-band photometry, with a library of stellar popu- 
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Fig. 3. The observed relationship between stellar mass and luminosity for galaxies in the lOK sample, within the three redshift 
ranges studied in this paper The left panel shows an aspect of the galaxy bimodality, with red galaxies more massive and brighter 
than blue ones. 



lation models based on iBruzual & Chariot! (l2003h . We used the 
code Hyperzmass, a modified vers i on of the photometric red- 
shift code Hyperz jBolzonella et al.Ll2000l) . The typical error on 
stellar masses is ~0.2 dex. The method and accuracy of these 
measurements are fully described in Bolzonella et al. (2009) and 
|PozzettietalJ(f2009l) . 

We have thus constructed a set of mass-selected samples, 
containing galaxies more massive than a given threshold. We 
choose the same redshifts ranges as used for the luminosity se- 
lected samples. The properties of the selected subsamples are 
summarized in Tableland represented in Fig.|2] 

4.3. Mass completeness 

Due to the flux-limited nature of surveys like zCOSMOS (Iab < 
22.5) the lowest-mass samples are affected to varying degree by 
incompleteness related to the scatter in the mass-luminosity rela- 
tion (Fig.|3]l. This introduces a bias against objects which would 
be massive enough to enter the mass selected samples, but too 
faint to fulfill the apparent magnitude limit of the survey. These 
missed high mass-to-light ratio galaxies will be those dominated 
by low-luminosity stars, i.e. the red and faint objects. Clearly, if 
this is not accounted for in some way, it would inevitably affect 
the estimated clustering properties, with respect to a truly com- 
plete, mass-selected sample (iMeneux et al.Ll2008l) . It is therefore 
necessary to understand in detail the effective completeness level 
in stellar mass of the samples that we have defined for our anal- 
ysis. 

iMeneux et al. I (I2008h have used 2 different methods to ex- 
plore and quantify the completeness limit in stellar mass as a 
function of redshift. The first is based on the observed scatter in 
the mass-luminosity relation, obtained from the data themselves 
and extrapolated to fainter fluxes. The second instead makes use 
of mock survey samples, under the hypothesis that they pro- 
vide a realistic description of the mass-luminosity relation and 
its scatter: the DLB40 set of mock survey catalogues that are 
complete in stellar mass are "observed" under the same condi- 
tions as the real data, i.e., selected at / < 22.5. The complete- 
ness is then simply defined, for a given redshift range and mass 




redshift z 

Fig. 4. Estimate of how the completeness in stellar mass changes 
as a function of redshift, due to the survey flux limit {Ij^b < 22.5). 
The shaded grey area and green contours describe the loci of 
constant completeness. They are derived from the DLB40 mock 
samples of 1 x 1 deg^ and defined as the fraction of observed 
{Iab < 22.5) galaxies over the total number in a given cell with 
size Az - 0.01 and log(M) > log(Mc,„). The red points super- 
imposed coiTespond to the actual data of the lOK sample. The 
yellow points and dotted line show the 95% M/L^atio com- 
pleteness level derived independently by Pozzetti et al.l (l2009l) 
(see text), directly from the observed data. The agreement be- 
tween the two estimates, from the data and from the Millennium 
mocks, is remarkable and adds confidence in the use of the sim- 
ulated samples. 



threshold, as the ratio of the number of galaxies brighter than the 
zCOSMOS flux limit over those at any flux. Interestingly, even 
if this method is model-dependent (in particular, on the prescrip- 
tion of galaxy formation used in the semi-analytical models). 
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Table 3. The completeness in stellar mass of mass-selected 
mock sub-samples reproducing the properties and selection cri- 
teria of our lOK data samples. Completeness is defined as the 
percentage of all galaxies above the mass limit which are actu- 
ally included in the sample. 



0.2 



redshift z 
0.5 0.8 



Sample 


Redshift 


Stellar mass {log{M/Me)) 


Completeness 




range 


range 




Ml.l 


0.2-0.5 


>9.0 


0.783 


M1.2 


0.2-0.5 


>9.5 


0.972 


Ml. 3 


0.2-0.5 


>10.0 


1.000 


M1.4 


0.2-0.5 


>10.5 


1.000 


M2.1 


0.5-0.8 


>9.0 


0.349 


M2.2 


0.5-0.8 


>9.5 


0.652 


M2.3 


0.5-0.8 


>10.0 


0.919 


M2.4 


0.5-0.8 


>10.5 


0.996 


M3.1 


0.8-1.0 


>10.0 


0.571 


M3.2 


0.8-1.0 


>10.5 


0.882 



this approach leads to similar completeness limits than the first 
one. The results of this second exercise are shown, as a func- 
tion of redshift and mass threshold and for a flux limit / < 22.5, 
in Fig. |4] Completeness is estimated in narrow redshift ranges 
(Az = 0.01) for diff'erent mass thresholds M^.,,, increasing from 
10** to 10"-^ h-^ Mq with a step of 10*'°' h'^ Mq. A large frac- 
tion of low-mass objects is clearly missed at high redshift. 

The yellow points and dashed line in Fig. |4] have been esti- 
mated from the observed scatter in the M/L relation of the data, 
and are defined at each redshift as the lower boun dary, Mmm(z) , 
inclu ding above it 95% of the mass distribution dPozzetti et al.L 
12009) . 

It is very encouraging to notice the very good agreement be- 
tween this independent estimation from the data and that based 
on the DLB40 set of mock catalogues. Table [3] summarises the 
completeness estimates derived from these mock catalogues for 
each of the 10 zCOSMOS galaxy samples defined in Table |2l 
The sample M2.1 shows the strongest incompleteness: 65.1% 
of the galaxies more massive than 10*^ h-^ Mq are fainter than 
/ = 22.5 at z=[0.5-0.8] and then, not included in our sample. In 
Sec. 16.31 we shall discuss the effects of this incompleteness on 
the galaxy clustering measurement. 



5. Estimating the two-point correlation function 

The two-point correlation function is the simplest estimator used 
to quantify galaxy clustering, being related to the second mo- 
ment of the galaxy distribution, i.e. its variance. In practice, it 
describes the excess probability ^(r) to observe a pair of galaxies 
at a given separation r, with respect to that of a random distribu- 
tion jPeebles, 1980). Here we shall estimate the redshift-space 
correlation function ^{rp,ji), which allows one to account and 
correct for the effect of peculiar motions on the pure Hubble 
recession velocity. In this case, galaxy separa tions are spht into 
the t angential and radial components, Vp and n dPavis & Peeblesl 
[1983; Fisher etal., 1994). 

The real-space correlation function ^s(r) can be recovered 
by projecting ^(r^, n) along the line-of-sight, as 



(2) 



For a power-law correlation function, ^^(r) = {rlr^Y^, this inte- 
gral can be solved analytically and fitted to the observed Wp{rp) 
to find the best-fitting values of the correlation length ro and 
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Fig. 5. Overall radial distribution of the zCOSMOS lOK sample, 
compared to three different smoothed distributions. These are 
obtained by filtering the observed data with a Gaussian kernel of 
increasing cr=150, 250 and 450 h"' Mpc. The first two smoothed 
curves retain information of the two large structures located at 
-1000 and -1800 h ' Mpc along the line-of-sight, while the 
third one overestimates the number density of galaxies at low 
and high redshift. 
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Fig. 6. The radial distribution of the luminosity selected sam- 
ple L2.1 is compared to a smoothed curved (with a kernel of 
cr-A50 h"' Mpc - dashed curve) and a radial distribution gen- 
erated from the integration of the luminosity function (solid 
curve). The latter is consistent with that expected for such galaxy 
sample. 
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slope J (e.g. jDavis & Peebleslll983h . In computing Wp{r^, a fi- 
nite upper integration limit has to be chosen in practice. Its value 
has to be large enough as to include most of the clustering signal 
dispersed along the line of sight by peculiar motion. However, it 
must not be too large, to avoid adding j ust noise, which is dom- 
inant above a certain n. Previous works jPoUo etal.L [20051) have 
shown that, for similar data, the best results are obtained with 
an integration limit n^ax between 20 and 40 h"' Mpc. Our tests 
show that the scatter in the recovered Wpixp) is obtained using the 
lowest value in this range. This can introduce a 5-10% under- 
estimate in the recovered large-scale amplitude, which can be 
accounted for when fitting a model to Wpivp). In the following, 
we shall in general use Timax - 20 h^'Mpc and show examples 
of how the amplitude is biased by this choice for the real data. 

To estimate ^{rp,n) from each galaxy sample, we use the 
standard estimator of Landv & Szalav ( 1993): 

^, , NR{NR-l)GG(rp,n) NR-lGR{rp,n) 

t(r„,n) H 1 (3) 

^ ' Nc(Nc-l) RR(rp,7T) Nc RR(rp,n) 

where Nq is the mean galaxy density (or, equivalently, the total 
number of objects) in the sample; A^^ is the mean density of a cat- 
alogue of random points distributed within the same survey vol- 
ume and with the same selection function as the data; GG{rp,ji) 
is the number of independent galaxy-galaxy pairs with separa- 
tion between rp and -i- drp and between ji and ji + dn; RR{rp, n) 
is the number of independent random-random pairs within the 
same interval of separations and GRijp, n) represents the num- 
ber of galaxy-random cross pairs. 

6. Observational biases and selection effects 

6. 1 . Correction of VIMOS angular foot-print and varying 
sampling 

To properly estimate the correlation function from the lOK 
zCOSMOS data, we need to correct for its spatial sampling rate, 
which is on average ~30%, but varies with the position on the 
sky due to the VIMOS foot-print and the superposition of mul- 
tiple passes (see Fig.lTJ. The correction s cheme used here is an 
evolution of that discussed in iPoUo et all ( 12005b . but with a sim- 
plified weighting scheme. The main differe nces of this sampl e 
with respect to the VVDS-Deep data used bv'P ollo et aLl(l2005l) . 
are that: (a) this sample is 1 .5 magnitudes brighter, (b) the spec- 
tra are taken with higher resolution, which produces longer spec- 
tra and thus less objects along the dispersion direction and, (c) 
there are as many as 8 repeated observations ("passes") cover- 
ing each point on the sky in the central area of the COSMOS 
field. The net result of these differences can be appreciated in 
Fig. [T] the sample is characterized by a well-sampled central re- 
gion, but also by rather sparsely sampled VIMOS pointings in 
the outskirts of the field. In particular, these external pointings 
clearly show target galaxies concentrated along rows. This ef- 
fect is produced by the significant length of the spectra on the 
CCD in medium-resolution mode: not more than 2 spectra can 
be aligned on top of each other on the detector in each quad- 
rant, which results into the observed two "stripes". This is sig- 
nificantly different from what happens i n the low-resolution ob - 
serving mode (as e.g. in VVDS-Deep, iLe Fevre et al.L l2005al) . 
where spectra are shorter and up to four of them can be packed 
along the same column on the CCD. 

We have tested three different algorithms to correct for the 
angular selection function of the survey, obtaining comparable 
results. Other weighting schemes use in particular the angular 



correlation functions of the lOK sample and the photometric 
catalogue to correct for the non uniform spatial sampling rate. 
Th ese methods are discus sed i n the parallel c lustering analyses 
bv lde la Torre et alj (|2009|) and lPorciani et al.l (|2009). In the lat- 
ter paper in particular, comparative tests of the three algorithms 
are presented. 

Since the sub-samples analyzed in this work are essentially 
volume-limited (above the luminosity/mass completeness lim- 
its), we do not need to apply any further minimum- varianc e 
weighting scheme (as e.g. the Jt, weighting. iFisher et aUll994b . 
This is normally necessary for purely flux-limited surveys in 
which the selection function varies significantly as a function of 
redshift, such that different parts of the volume are sampled by 
galaxies with different luminosities and number densities (e.g., 
lLietal.Ll2006h . 

6.2. Construction of reference random samples 

A significant source of uncertainty that we encountered in esti- 
mating two-point functions from our lOK sub-samples is related 
to the construction of the random sample and in particular to its 
redshift distribution. We soon realized that the strongly clustered 
nature of the COSMOS field along the line of sight, with sev- 
eral dominating structures at different redshifts, required some 
particular care as not to generate systematic biases in the ran- 
dom sample. These superclusters are already evident as vertical 
stripes in Fig. |2] and even more clearly in the redshift histogram 
of Fig.|5] Note the big "walls" at z = 0.35, 0.75 and 0.9, which 
are also c learly identified by the density field reconstruction of 
iKovac et al. (2009). 

A standard way to generate a random redshift coordinate 
accounting for the radial selection function of the data uses a 
Gaussian-filtered version of the data themselves. This is nor- 
mally obtained using smoothing kernels with a dispersion cr (in 
co-moving coordinates) in the range 150-250 h"' Mpc. The 
results of applying this technique to the current lOK data are 
shown in Fig.|5] One notes how for smoothing scales of 150 and 
250 h"' Mpc the curves still retain memory of the two largest 
galaxy fluctuations. These are erased only when a very strong 
smoothing filter (450 h ' Mpc) is adopted. However, in this case 
the smoothed curve is unable to follow correctly the global shape 
of the distribution, over-estimating the number density in the 
lowest and highest redshift ranges. The situation for our specific 
analysis, however, is somewhat simpler than this general case. 
Our luminosity-limited or mass-limited samples are in princi- 
ple "volume-limited", i.e. samples that - if properly selected - 
should have a constant density within the specific redshift bin. 
One such case is shown in the zoom of Fig.|6] where the redshift 
distribution in the range z = [0.5, 0.8] is plotted. 

An alternative way to generate the radial distribution of the 
random sample is to integrate the galaxy luminosity function 
(LF) in steps along the redshift direction, computing at each 
step a value for the density of galaxies expected at that redshifts. 
Ideally, the LF can be measured from the sample itself and would 
include any detected evolution of its parameters. This is what we 
did here, using the evolving LF parameters presented in the com- 
panion dedicated paper (Zucca et al. 2009). The red dashed line 
in Fig. |6] shows the result one obtains if smoothing with a kernel 
of 450 h"^ Mpc (dashed), compared to that obtained from the 
integration of the LF (solid). The latter is fully consistent with 
that expected from a truly volume-limited sample with the given 
selection criteria, with the number of objects increasing as the 
square of the radial comoving distance. 
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Fig. 7. The effect of stellar mass incompleteness on the measured 
Wp{rp), estimated from the 1x1 deg- Millennium mock samples. 
The figure shows the average of the quantity R over 40 mock 
samples as a function of rp. R is defined as the ratio of the esti- 
mates of Wp{rp) with and without the flux cut at Iab = 22.5, i.e. 
for a sample mimicking the lOK selection and a sample 100% 
complete in stellar mass. 



6.3. Effect of mass incompleteness on Wp(rp) 

As discussed in Sec. 14.31 when we constructed our mass-limited 
samples, a fraction of galaxies more massive than the formal 
mass threshold are in fact lost due to the limiting Iab < 22.5 
flux cut of the survey. This becomes more and more important 
with increasing redshift. As we said this population of missing 
galaxies is inevitably dominated by red objects with high mass- 
to-light ratio (Meneux et al.. 2008,) . which are known to clus- 
ter m ore strongly than the average pop ulation (iMeneux et al.L 
l2006HCoil et all 120081: iMcCracken et al.Ll2008h . Having defined 
our clustering tools, we can now further extend the analysis of 
Sec. 14.31 and use the DLB40 mock samples to quantify directly 
the effect this has on the measured Wp(rp). We thus computed 
the statistics for each of the mocks, which are complete down 
to very small masses (~ 10^ Mq), with and without applying 
the apparent-magnitude cut. Clearly, we are making here a very 
strong hypothesis, i.e. that the simulated samples have intrinsic 
clustering properties (and their relation to the galaxy's M/L ra- 
tio), that are similar to those of real data. 

The ratio of these two estimates ("true" over "observed") av- 
eraged over the 40 mock catalogues is shown on Fig. |2l For a 
mass selection which is 100% complete within the given red- 
shift bin we would measure {R} = 1 at all separations. We can 
see that the only mass range for which this is strictly happening 
at any redshift is that with log(M/M0) >10.5. For smaller mass 
samples we see a clear reduction of the clustering amplitude. 
However, we can also see that for most samples the shape of 
Wp{rp) is distorted mainly only below < 1 h"' Mpc. Above this 
scale, the mass incompleteness introduces an amplitude reduc- 
tion up to ~ 20% in the worst cases. This will have to be consid- 
ered when comparing our measurements with models (although 
keeping in mind that these estimates come from simulated data, 
not from real observations). For general comparisons, however, 
the amount of amplitude reduction of Wp(rp) is typically negligi- 
ble on scales larger than ~1 h"' Mpc, given the statistical errors 
of the data measurements. 
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Fig. 8. Ratio of the diagonal errors on Wp(rp), obtained through 
the bootstrap resampling method to the "true" ones obtained 
from the variance of 24 mock catalogues. Filled and open sym- 
boles correspond to two different bootstrapping techniques, re- 
sampling respectively sub-volumes of the survey or single galax- 
ies. The former technique clearly provides a standard deviation 
which is closer to the "true" one obtained from repeated mea- 
surements. 



7. Systematic and statistical errors on correlation 
estimates 

The derivation of realistic errors on the galaxy correlation func- 
tion has been the subject of debate since its early measurements 
(see e.g.. Bernstein. JL994). In particular, it is well-known that the 
measured values of the two-point correlation function on differ- 
ent scales are not independent. This means that, e.g., the bins of 
Wp(rp) have a degree of correlation among them, that needs to be 
taken into account when fitting a model to the observed values. 
This can be done if we are able to reconst ruct the N x N cova ri- 
ance (or correlation) matrix of the N bins (Fisher et al.LlT994b . 

In a recent paper, iNorberg et al.l ( 12008 ^ compare in detail 
three different methods for estimating the covariance matrix of a 
given set of measurements. These use: (a) the ensemble variance 
from a set of mock catalogues reproducing as accurately as pos- 
sible the clustering properties and selection function of the real 
data: (b) a set of bootstrap resamplings of the volume contain- 
ing the data, and (c) a so-called jack-knife sub-set of volumes of 
the survey. In this latter case, the survey volume is divided into 
Nv sub-volumes and the statistics under study is re-computed 
each time excluding one of the sub-par ts. In the "block-wise" 
incarnation of the bootstrap technique (P orciani & Giavaliscol 
2002, method "b"), instead, sub-volumes are selected each 
time with repetition, i.e. excluding some of them, but counting 
two or more times some others as to always get a global sample 
with the same total volume. We note however that historically 
there are two possible ways to resample interna lly the data set. 
The classical "old" bootstrap (iLing et al.L 1 19861) entailed boot- 
strapping the sample "galaxy-by-galaxy". This means picking 
randomly each time a sample of Nc galaxies among our data 
set of Nc galaxies, allowing repetitions. In this way, within one 
bootstrap realization a galaxy can be selected more than once, 
while some others are never selected. This technique has been 
shown to le ad in genera l to so me under-estimation of the diag- 
onal errors dFisher et al.L 119941) . Here we shall test directly also 
this aspect. 

The advantage of using mock samples is that, under the as- 
sumption that these are a realistic realization of the real data. 
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they allow us to obtain a true ensemble average and stan- 
dard deviation from samples with the same size as the data 
sample, including both Poissonian noise and cosmic variance. 
Unfortunately, the covariance properties derived from mock 
samples are not necessarily a good description of those of the 
real data, thus making the use of the derived covariance matrix 
(e.g. in model fitting) doubtful. Conversely, depending on the 
sample size, jack-knife or volume-bootstrap covariance matrices 
can exacerbate peculiarities of some sub-regions, again not rep- 
resenting adequately the true covariance properties of the data. 

For the present investigation, we spent considerable effort 
to understand how to best estimate a sensible covariance ma- 
trix for our Wp{rp) measurements. The available mock samples 
were crucial as to allow us to perform direct comparisons of 
the performances of the different techniques. After some ini- 
tial attempts, we excluded the jack-knife method because of the 
limited size of the survey volume. We then performed a direct 
comparison of the covariance matrices derived through the boot- 
strap technique and from the KW24 mock catalogues. For the 
bootstrap method, we decided to test directly how galaxy- and 
volume-bootstrap were performing. We concentrated on the red- 
shift range z= [0.5-0. 8] selecting simulated galaxies brighter than 
Mb - 5log(h) < -19.5 - z. After computing the correlation func- 
tion Wp(rp) for all 24 mock samples, we constructed for each 
of them: (a) 100 galaxy-galaxy bootstrap samples and (b) 100 
volume-volume bootstrap samples. In the latter case, we con- 
sidered 8 equal sub-volumes, defined as redshift slices within 
the redshift range considered. The number of sub-volumes was 
choosen as the best compromise between having enough of them 
and not having too small volumes. With this choice, their volume 
is ~ 1.4 X 10^ h"^ Mpc^ for the samples with z=[0.5-0.8] and 
z=[0.8-1.0] and ~ 0.6 x 10^ h"^ Mpc^ for z=[0.2-0.5]. The two 
bootstrap techniques lead to a total of 4800 samples and corre- 
sponding estimates of Wp{rp). We then calculated the covariance 
(and correlation) matrices for each of these two cases along with 
the one derived from the correlation function of the 24 mocks 
themselves. 

In Fig. [8] we show a comparison of the standard deviations 
derived from the two bootstrap techniques, to that derived from 
the 24 mocks. In each case, these values correspond by defini- 
tion to the square root of the diagonal elements of the covariance 
matrix. In the plot we show the mean (over the 24 mocks) of the 
ratio of cr„,^^ from the bootstrap to the "true" one from the ensem- 
ble of mock surveys. This shows clearly how the r.m.s. values 
obtained with the single-galaxy bootstrap grossly underestimate 
the true variance, up to one order of magnitude on large scales. 
Bootstrapping by volumes produces a better result, providing a 
realistic estimate of cr„,^ between 0.1 and 1 h"' Mpc, and a 20- 
25% under-estimate on larger scales. 

Each element of the the correlation matrix r,j is obtained 
from the corresponding element of the covariance matrix crij 
as rij = cTijl ^a-iiO-jj. By definition, the off'-diagonal terms of 
the coiTelation matrix will then range between -1 and 1, indi- 
cating the degree of correlation between diff'erent scales of the 
function Wpivp). Considering the redshift range z=[0.5-0.8], we 
show in Fig. |5]the mean of the 24 correlation matrices derived 
resampling the galaxies (left panel), or resampling 8 equal sub- 
volumes (center), 100 times each. These are compared to the cor- 
relation matrix directly derived from the 24 mocks (right panel). 
The first case shows a mainly diagonal coiTelation matrix where 
off'-diagonal terms are mostly noise. In the second case they in- 
stead decrease smoothly from 1 to as a function of bin separa- 
tion. The matrix derived from the 24 mocks shows high correla- 
tion at all scales. 



Table 4. The 5 main eigenvalues of the correlation matrix de- 
rived with the bootstrap resampling, respectively of galaxies 
(first column) and sub-volumes (col. 2), and from the ensemble 
variance of the 24 mocks (col. 3). For the two latter cases, each 
mock sample is used in turn as "data" and the reported eigenval- 
ues are the obtained as the average over the 24 mocks. 



eigenvalue 


bootstrap 


mocks 




galaxies 


volumes 


catalogues 




3.87784 


8.28193 


11.82062 


^2 


1.95834 


2.17866 


0.15708 


^3 


1.36841 


0.71635 


0.01794 


^4 


1.11316 


0.38252 


0.00436 


^5 


0.91683 


0.20340 


0.00000 



In order to compare directly the properties of the correla- 
tion matrices derived with the 3 methods, we diagonalize all of 
the 24 -H 24 H- 1 matrices by computing their principal compo- 
nents and the amplitudes of the corresponding eigenvalues /I,- 
(/ = 1 - 12). Note that the sum of the eigenvalues of a correlation 
matrix is always equal to its dimension, i.e. 12 in our case. We 
report in Table |4]the values of the five main eigenvalues obtained 
with the 24 mocks (first column) compared to the averages over 
the 24 mocks of those obtained with the 2 resampling methods. 
The numbers show that the correlation matrix derived from the 
24 mocks contains essentially four principal components and is 
mostly dominated by one of them. This indicates a strong cor- 
relation in the data. The bootstrap matrices, instead, show more 
than 5 non-negligible components, with the fifth one being of 
the same order of magnitude of the second in the mock matrix. 
This implies a lower correlation. We note, however, that volume 
resampling tends to produce a matrix whose structure is closer 
to that of the mocks, with 1-2 dominant components. This is an- 
other indication of how volume bootstrapping, although also not 
reproducing perfectly the intrinsic covariance properties of the 
sample, provides a better estimate of the variance and correla- 
tion in the data with respect to galaxy-galaxy bootstrap. 

These experiments are extended and further d iscussed in our 
parallel accompanying papers, in particular by Porcia ni et al.l 
(2009). The bottom-line result of our extensive investigations is 
that volume bootstrap, if enough resamplings are used, provides 
a sufficiently good reconstruction of the intrinsic covariance ma- 
trix of the data set. This is obtained at the expenses of a slightly 
less accurate account of cosmic variance on large scale, with 
respect to what can be obtained from the scatter among mock 
samples, where wavelengths longer than the survey size can be 
sampled. However, we have shown (Fig. [SJ that this effect on 
scales ~ 10 h"' Mpc is limited to ~ 20%. 

8. Dependence of galaxy clustering on luminosity 

8.1. Luminosity dependence at fixed redshift 

Figure [TO] shows the projected correlation function Wpijp) es- 
timated for our 9 luminosity selected sub-samples at different 
redshifts. Error bars correspond to the Icr dispersion provided 
by 200 volume bootstrap resamplings, as extensively discussed 
in Sect. |7] 

We note that no clear dependence on luminosity is observed 
within any of the three redshift ranges. Also, in the shape of 
Wp{rp) there is some hint for the usual "shoulder", i.e. a change 
of slope around 1 h ' Mpc, but no clear separation between the 
expected 1-halo term on small scales and the 2-halo component 
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Fig. 9. Mean of the 24 correlation matrix derived resampling the galaxies of eack KW24 mocks {left), or resampling 8 equal sub- 
volumes {center). These are compared to the correlation matrix derived directly from the KW24 mocks {right). The redshift range 
considered here is z=[0.5-0.8]. The averaging over the 24 realisations of the 2 left matrices suppress the negative off-diagonal terms 
which are sometimes present for a given mock catalogues. Correlation coefficients are then color-coded from to 1. 
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Fig. 10. Projected correlation function Wp{rp) measured as a function of galaxy luminosity within three redshift ranges. No signifi- 
cant dependence on luminosity is observed within the explored ranges. Note the very flat slope of Wp{rp) in the central redshift bin, 
compared to the two other slices. 



above this scale (see the Introduction for definitions). In particu- 
lar in the intermediate redshift bin, all sub-samples show a rather 
flat large-scale slope, with no evidence for the usual breakdown 
above ~ 2 h"' Mpc. 

To try and understand the origin of the observed flat shape, 
it is interesting to look directly at the contour plots of the bi- 
dimensional redshift-space correlation function ^{rp,ji). These 
are shown in Fig. [TT| for the three luminosity-selected sub- 
samples LI .4, L2.2 and L3. 1 (see Table[T]for definitions), that in- 



clude galaxies brighter than Mg - 5log{h) < - 19.5 - z. The three 
contour plots show some interesting features. First, one clearly 
notices the much stronger distortion along the line of sight tt, in 
the central panel. At the same time, in the same redshift range 
a much more extended signal is also observed along the per- 
pendicular direction rp. It is tempting to interpret both these ef- 
fects as produced in some way by the two dominating structures 
that we evidenced in Fig.|6] The excess signal along the line of 
sight is very plausibly due to the distortions by "Fingers of God", 
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Fig. 11. Iso-correlation contours of ^{rp, n) (here smoothed with a gaussian kernel) for galaxies brighter than Mb - 5log(h) < 
-19.5 -z, computed in 3 redshift ranges. The amplitude is color-coded according to the side bar, while the black contour corresponds 
to ^{rp, n) - 1. White values correspond to ^{rp, n) < 0.4. 



due to an anomalous number of virialized systems (groups and 
clusters) within these structures. At the same time, the extension 
along rp is indicating that there is also an excess of pairs per- 
pendicularly to the line of si ght, with respect to an isotropic dis- 
tribution. In fact, we know jScoville et al.L l2007bl: iGuzzo et al.L 
[2007) that the large-scale structure at z ^ 0.73 extends over a 
large part of the COSMOS area. This evidently biases the ob- 
served number of pairs along rp, for simple geometrical reasons. 
We cannot exclude that part of the large-scale compression ob- 
served in ^(fp, n) is also generated by an excess of galaxy infall 
onto t his structure, th us producing what is known as the Kaiser 
effect (iKaiseil IT981 . This effect is proportional to the growth 
of structure (see e.g.. lGuzzo et al.L 1200 8, for a recent direct esti- 
mate at similar redshift) and can be extracted when the underly- 
ing clustering can be assumed to be isotropic. In this case it is in 
practice impossible to disentangle this dynamical distortion from 
the geometrical anisotropy generated by having one dominating 
structure elongated perpendicularly to the line of sight. 

The flatter shape in w„( rp) in Fig.[TOl is also consistent with 
the overdense samples of Abbas & Shethi (|2007|) . who noticed 
not only a higher amplitude for the most overdense (10% and 
30%) samples of mock and SDSS galaxies, but also a flattening 
in the correlation function with respect to the full sample. This 
is another line of evidence in favour of the hypothesis that the 
zCosmos field is centered on an overdensity. 

The plots of Fig. [TT| also show explicitly the reasons for our 
choice of nmax - 20 h"' Mpc as the upper integration limit 
in the computation of Wp{rp): this value provides a reasonable 
compromise between including most of the signal and exclud- 
ing the noisiest regions in the upper part the diagrams. In the 
central redshift bin, however, some real clustering power may 
still be present above this scale, for small r^'s. In Fig. [T2l we 
show directly how Wp{rp) changes, when n^ax is extended from 
20 to 30 h"' Mpc. We see that, somewhat counter-intuitively, 
below 1 h ' Mpc no extra amplitude is gained, while - as indi- 
cated by the mock experiments (see Sect. |5]) - the scatter is in- 
creased. Conversely, one can see the slight scale-dependent bias 
on the amplitude at larger separations, which gets up to ~ 10% 
at 15 h"' Mpc when increasing n,„ax- 
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Fig. 12. Sensitivity of the projected function Wp{rp) to the upper 
integration limit n,„ax, for one luminosity-selected sample in the 
central redshift bin. 



8.2. Redshift evolution at fixed (evolving) luminosity 



In section 14.11 we have discussed how our luminosity selection 
has been devised as to account for the average evolution in the 
luminosity of galaxies, assuming this to be the dominant ef- 
fect in modifying the mean density of objects at a given lu- 
minosity. Under this assumption, it is then interesting to test 
how Wp(rp) for galaxies within the same (de-evolved) luminos- 
ity interval changes with redshift. This is shown in Fig. [13] for 
Mb - 5log(h) < -19.5 - z, i.e. for the same three samples for 
which ^{rp,n) is plotted in Fig. [TT] No coherent evolution of 
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Fig. 13. Evolution of the projected function Wp(rp) of galaxies 
with Mb - 5log(h) < -19.5 - z between redshift z - 0.2 and 
z= 1. 



the amplitude and shape of the projected correlation function 
is observed among the three samples, characterized by mean 
redshifts 0.37, 0.61 and 0.91. The three curves are consistent 
with each other within the error bars, with - as expected - the 
intermediate-redshift bin (triangles) showing a systematically 
higher amplitude than the other two. Again, this is easily in- 
tepreted as a local effect, resulting from the extreme large-scale 
clustering observed in this redshift range. On the other hand, 
the overall lack of apparent amplitude evolution of luminous 
galaxies with red shift is consistent with prev ious results from 
the VVDS- Deep dPollo et al .'. '2006'). DEEP2 ( Coilet al.Ll2006h 
and SDSS jZehavi et all 2005) surveys for galaxies brighter 
than ~ M*. The only evolutionary effect evidenced in particu- 
lar by the VVDS data is a steepening of the small-scale slope of 
H'p(yp) (i.e. t he 1-halo term) for high-redsh ift luminous galaxies 
dPollo et al.Ll2006l) . lMcCracken et alj 12008 ) also show alack of 
clustering amplitude evolution for a large sample of luminous 
(-22 < Mb - 5log{h) < -19) galaxies based on accurate photo- 
metric redshifts in the CFHTLS Deep fields; interestingly, they 
show that such invariance is maintained also when splitting the 
sample into early- and late-type galaxies. 



9. Dependence of galaxy clustering on stellar mass 

The relation of clustering properties to galaxy stellar masses is 
in principle more informative and straightforward to interpret 
as stellar mass is a more fundamental physical parameter than 
luminosity. 



9.1. Mass dependence at fixed redshift 

Also in the case of stellar mass dependence, it is interesting to 
look at the shape of the full correlation function ^(rp,7i) in the 
three redshift ranges. We show in Fig. [141 the result obtained 
for the 3 samples Ml. 3, M2.3, M3.2 (see Table|2|, that include 



galaxies more massive than 10^° h^^ Mq. Also in this case the 
central panel is significantly different from the other two, with 
^{rp, n) remaining positive out to much larger scales in both rp 
and 71 directions. We note that the small-scale stretching along 
71 seems to be less extended than that obtained for the corre- 
sponding luminosity-selected sample, although it is hard to say 
whether this difference is significative. 

Figure [15] shows the projected correlation function Wp{rp) 
of the 10 mass-selected samples. The plotted points are not cor- 
rected for the residual stellar mass incompleteness (see Sec. |6.3t . 
Errors are estimated as in the luminosity case using 200 boot- 
strap resamplings of 8 equal sub-volumes of each dataset. The 
figure shows a weak mass dependence of clustering in the low- 
and high-redshift bins, in particular at small separations. At the 
same time, a strong dependence at all separations is evident in 
the intermediate redshift slice. In this range the slope of Wp(rp) 
remains extremely flat out to the largest explored scales, even 
more strongly than in the luminosity-selected cases. Finally, we 
note also that in the low- and high-redshift bins there is evidence 
for a steeper "1-halo term" contribution at < 1 h"' Mpc (with 
no clear indication for an evolution in redshift of the transition 
scale to the 2-halo term). Conversely, the central redshift range 
seems to be characterized by the same, flat power-law shape 
down to 0.2 h ' Mpc, where a sudden steepening is then ob- 
served. The slope below 0.2 h ' Mpc seems to depend directly 
on the limiting mass, with more massive galaxies showing a 
steeper correlation function. In summary, no clear overall trend 
can be evidenced among the three redshift ranges, with the cen- 
tral volume again displaying peculiar clustering properties that 
apparently dominate over any possible cosmological effect. 

9.2. Ciustering evoiution at fixed stellar mass 

It is also interesting to compare directly the evolution of Wp(rp) 
with redshift, when a specific class of stellar mass is selected. As 
mentioned earlier in this section, this is particularly interesting, 
as in principle it does not require accounting for strong galaxy 
evolutionary trends as in the case of luminosity. We are assum- 
ing here that stellar mass does not significantly grow between 
z ~ I and z ~ 0.2, which we know is only partially true. A fac- 
tor of ~ 2 growth in stellar mass is in fact expected on average 
between z - I and z = 0, which however would have little effect 
on the estimated correlation function, if taken into account. In 
Fig. [16] we show Wp{rp) computed for the same three samples 
with logiMIMo) > 10 of Fig. [14] Similarly to the luminosity 
samples, we do not see any clear evolution with redshift of the 
amplitude and shape of the projected correlation function. The 
three curves are consistent with each other within Icr. 



10. Comparison with independent measurements 
and models 

10.1. Redshift evolution of Wp(rp) 

An accurate z ~ reference measurement of Wp( rp) as a functio n 
of stellar mass has been obtained by the SDSS (iLi et al.Ll2006l) . 
Meneux et al. (2008) do find evidence for evolution of the am- 
plitude of Wp{rp) for galaxies less massive than 10'°'^ h^^ Mq, 
when comparing this to the measurements from VVDS-Deep 
at z ~ 0.9. The SDSS and zCOSMOS ste llar masses wer e de- 
rived with the same initial mass function ( IChabrieJ, 12003 *) and 
normalised to h=l. They are directly comparable. The SDSS 
clustering meas urements were obtained within differential stel- 
lar mass ranges (iLi et al.Ll2006l) . while ours correspond to galax- 
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Fig. 14. Example of full redshift-space correlation function ^{rp, n) (here smoothed with a gaussian kernel) for the sub-sample 
containing galaxies more massive than 10'" h^^ Mq computed in the same 3 redshift ranges as in Fig.[TT] ^{rp,n) is computed in 
cells of 1 h ' Mpc side in both rp and n and the iso-correlation levels are color coded according to the side bar. The thick black 
contour corresponds to ^(r^, n) - \ and the white values to ^(fp, jf) < 0.4. Note how the central panel (z = [0.5 - 0.8]) shows extra 
power in both directions, perpendicular and parallel to the line-of-sight. The small-scale "Finger-of-God" effect is however less 
pronounced than in the case of the luminosity-selected sample of Fig.fTTI 
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Fig. 15. The projected correlation function Wpivp) as a function of galaxy stellar masses in the zCOSMOS lOK sample, within three 
redshift ranges. 



ies more massive than a given threshold. However, from Fig. |2] 
we see that the zCOSMOS sample includes a very small number 
of galaxies more massives than lO" Yc^ Mq, due to the much 
smaller volume when compared to the SDSS. Any of our mass- 
selected samples has therefore, in practice, an upper bound at 
this value of mass. This implies that we can coherently com- 
pare two of the SDSS measurements of Wp(rp) (for their galaxy 
samples with stellar masses in the ranges [10.0 - 10.5] and 



[10.5 - 11.0]) to those from our samples M3.1 and M3.2, that 
include galaxies more massive than 10'" and lO'"'^ h"^ Mq re- 
spectively, within the redshift range z= [0.8- 1.0]. This compari- 
son (Fig. [TtT i does not show a clear evolution with redshift. For 
both samples, the large-scale amplitude of Wp{rp) is virtually the 
same as in the local SDSS samples. Considering simple evolu- 
tion of structures, this implies that the bias for galaxies more 
massive than 10'" h"^ Mq, has evolved significantly between 
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Fig. 16. Evolution of the projected function Wp(rp) of galaxies 
with logiMIMo) > 10 between redshift z - 0.2 and z-l. 



z ~ \ and today, as to keep their apparent clustering ampli- 
tude substantially unchanged. This implies in practice that the 
bias b{z) must evolve in a way such that b{Qi)D{Qi) ^ b{z)D(z), 
where D{z) is the linear growth rate of density fluctuations. In the 
standard model, this implies that at the approximate mean red- 
shifts of our redshift bins, z - 0.35, 0.75, 0.9, the bias of massive 
galaxies must have been res pectively 1.2, 1 .44 an d 1.53 times 
its value at the current epoch jMeneux et al] (l2008l) observed the 
same effect at < z >~ 0.8 in the VVDS data but only for galax- 
ies more massive than ~ 10'"^ h"^ Mq, with lower-mass objects 
showing a weaker bias evolution. A non evolution of the cluster- 
in g of the mo st massive galaxies was also noticed in the NDWFS 
jBrown et all 120081) and 2SLAQ surveys dWake et al.Ll2008l) . 



10.2. Observed and predicted shape ofwp(rp) af 0.5 < z < 1 

The only available measurement of clustering as a function 
of stellar mass at redshifts comparable to thos e explored by 
our s ample is that from the VVDS-Deep survey (M eneux et al.L 
l2008i) at 0.5 < z < 1.2. VVDS-Deep goes 1.5 magnitude deeper 
(although over a smaller area of ~ 0.5 deg^), which allows the 
analysis to be extended beyond z - 1 . To provide a qualitative, 
yet meaningful comparison of these two data sets, we can re- 
compute the correlation function for the lOK data within the 
largest usable redshift range overlapping with the VVDS inter- 
val, i.e. [0.5-1.0]. We applied the same stellar mass selection lim- 
its, keeping in mind the residual incompleteness that will affect 
the highest redshift part of the sample. The result is shown in 
Fig. [13 where the VVDS and zCOSMOS mass-selected samples 
are directly compared. The difference in shape and amplitude in 
the Wp(rp) derived from the two data sets is rather striking. The 
zCOSMOS points show in general a much flatter relation than 
those from the VVDS. The amplitudes for a given mass selec- 
tion also seem to be incompatible at several standard deviations 
between the two samples, especially above 1 h ' Mpc. 
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Fig. 18. Direct comparison of the dependence of clustering on 
stellar mass in the VVDS-Deep and zCOSMOS samples, over 
a similar redshift range. This plot evidences the significant in- 
trinsic difference between the two surveys, with the zCOSMOS 
sample showing a significantly larger clustering amplitude for all 
stellar masses. In addition, the much flatter shape of Wp(rp) indi- 
cates the predominance of coherent structure perpendicularly to 
the line of sight, due to the "walls" discussed in the text. 



10.3. Comparison to analytic and semi-analytic models 

It is at this point relevant to compare the available observations 
from zCOSMOS and VVDS with model predictions in a stan- 
dard ACDM scenario. We ca n do this in two way s. We first used 
the HALOFIT public code (ISmifli et all 120031) . fliat uses the 
halo model to compute the expected non-linearly evolved power 
spectrum at z = 0.8, that we take as a reasonable mean redshift 
for the two samples (our conclusions would not differ at all if 
predictions for z = 0.7 or 0.9 were used). The corresponding pro- 
jected function Wp{rp) is then computed by Fourier-transforming 
the power spectrum and projecting the resulting real-space cor- 
relation function. The result gives the expected Wp(rp) of the 
mass density field at z = 0.8. Secondly, we can compute the 
expectation value and the scatter expected in the same redshift 
range (0.5 < z < 1) for Wp(rp) using the available semi-analytic 
mock surveys built from the Millennium run. To this end, we 
use the DLB40 mocks for which we have full control over stel- 
lar masses, selecting simulated galaxies with I^b ^ 22.5 and 
log(M/M0) > 10, reproducing the sampling rate of the lOK 
data. In Fig. [19] we plot both the HALOF/T prediction for the 
dark matter (lower solid line) and for an arbitrarily biased pop- 
ulation of halos with b - V2!6 ~ 1.61 (dashed line), together 
with the BDL40 mean Wp(rp) (blue, lighter solid line) and the 
corresponding Icr and 3cr scatter corridor from the 40 mock sur- 
veys (solid- and dotted-dashed areas, respectively). We note as 
a consistency check the rather good agreement between the ana- 
lytic HALOFIT result and the expectation value from the full n- 
body plus semi-analytic simulation. On these model predictions, 
we overplot the corresponding zCOSMOS and VVDS estimates. 
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Fig. 17. The measured projected function Wp(rp) for galaxies at z=[0.8-1.0] with logiM/Mo) ~ [10 - 11] (left) an d logjMIM^:,) ~ 
[10.5 - 11] (right) from the zCOSMOS survey (filled diamonds), compared to the z ~ 0.1 estimate by the SDSS jLi et all l2006l) 
(blue curve). 



The zCOSMOS points agree rather well with the models (at bet- 
ter than 68% confidence) in both shape and amplitude on scales 
smaller than 1 h ' Mpc. On larger scales, however, the observed 
Wpifp) would require a strongly scale-dependent bias to be com- 
patible with the model predictions. Such a scale-dependence 
would also behave oppositely to what models and very general 
considerations suggest, implying a bias which grows with scale, 
rather than declines. From the plot we see in fact that the lOK 
data are compatible with b - 1 .6 on small scales, but would re- 
quire b ^ 2.45 on 10 h"' Mpc scales. The shaded area shows 
that this large-scale excess is marginally compatible with the 
model predictions, representing a very strong positive fluctua- 
tion. A few percent of volumes this size would show such a high 
clustering amplitude (on these scales and for this kind of galax- 
ies), in a ACDM Universe. 

It is interesting to note, at the same time, how the VVDS 
measurements lie on the opposite side of the distribution, at 
about 1.5-2crfrom the mean, but with a shape which is compati- 
ble with the model prediction over the whole ran ge (coiTespond- 
ing to a linear bias b ~ 1.2). Comparison with lAbbas & ShethI 
(120071) also suggests that the zCOSMOS field is centered on an 
overdensity, whereas the VVDS field is a less significant under- 
density. 

These results show how a full HOD model fitting to the 
Wp{rp) measured from the lOK data - that we originally planned 
to include in this paper - would add no meaningful information 
to the cuiTent analysis. Our first experiments with HOD mod- 
els based on the universal halo mass function indicate that rather 
unrealistic sets of parameters are required to reproduce the ob- 
served function. An interesting possibility would be to use in 
such modelin g a halo mass ftinctiori that d epends on local envi- 
ronment (e.g. lAbbas & ShethI |2005LI2006|) . to take into account 



the evidence that a large part of this sample is dominated by an 
overdensity. We plan to explore this possibility using the larger 
20K zCOSMOS sample which is now nearly complete. 

11. Discussion 

Together w ith previous analyses d McCracken et al.L 120071 : 
iKovac et"an, [2009), these results suggest that a significant frac- 
tion of the volume of Universe bounded by the COSMOS field 
is indeed characterized by particularly extreme density fluctu- 
ations. We have seen how, in statistical terms, these seem to 
lie at the 3cr limit of the distribution of amplitudes expected 
in volumes of a few 10^ h"^ Mpc^. We should consider, how- 
ever, that these conclusions are drawn from measurements that 
are strongly affected by the angular distribution of structure. The 
McCrac ken et alj (120070 result is based on the angular correla- 
tion function, while here we have studied the projected function 
Wp(rp). The latter, although making use of the redshift infor- 
mation is in practice a clustering measure dominated by galaxy 
pairs lying almost perpendicularly to the line of sight. The under- 
lying assumption when measuring Wp(rp) is that the geometri- 
cal distribution of structures within the sample being analyzed is 
completely isotropic. In other words, that there are superclusters 
aligned along several directions, such that the only remaining 
radial signal is produced by galaxy peculiar velocities. The very 
reason of using Wp{rp) is indeed to get rid of the distortions in- 
troduced in the shape of ^(s) (the redshift-space, angle-averaged 
coiTelation function) by galaxy motions. If this is true, and only 
in this case, then Wp(rp) is fully equivalent to an integral over 
the real-space correlation function ^(r) and therefore carries the 
same cosmological information. However, if, as in the case we 
have encountered here, there is one or more dominating struc- 
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Fig. 19. The zCOSMOS (solid circles) and VVDS-Deep (open 
circles) Wp{rp) of galaxies with mass larger than 10'" /i"^ Mq and 
< z >~ 0.8, compared to model predictions. These include the 
non-linear mass projected c orrelation function computed using 
HALOFIT (lower solid line, ISmith et al.L l2003h . and, for refer- 
ence, the corresponding Wpifp) for a population of halos with 
bias - 2.6 (dashed curve). The latter curve is a very good de- 
scription of the full non-linear Wp(rp) (light solid blue line), ob- 
tained averaging the forty DLB40 mocks from the Millennium 
run after applying the same sampling, magnitude and mass se- 
lections of the lOK sample. The shaded areas (thick and thin 
shaded) give, respectively, the Icr and 3cr confidence corridors 
around the mean. The large-scale zCOSMOS Wp{rp) in this red- 
shift range is thus marginally compatible with a rare, ~ 3cr pos- 
itive fluctuation in a standard ACDM Universe. 



tures extending preferentially along one direction, then the use 
of Wp(rp) to infer cosmological information is inappropriate. 

One may thus wonder whether more robust cosmological in- 
formation could instead be inferred looking directly at the sim- 
plest, angle-averaged redshift-space correlation function ^{s). 
The expectation is that the average over all directions reduces 
the weight of the excess pair counts produced by just a few struc- 
tures oriented along one preferential direction. In such case any 
analytic modelling (e.g. with HOD models) should also include 
an appropriate model for the linear and non-linear redshift distor- 
tions (IScoccimarro, 2004; Tinker et al., 2007). More simply, we 
can use the available mock samples in redshift space to compute 
the non-linear redshift-space ^{s) and its variance and compare 
it to the data, as we did for Wp(rp). In Fig.|20]we first plot f(s) 
for the lOK sample, computed for the usual four mass ranges in 
the broad redshift range 0.5 < z < 1 . The four data sets show 
a smooth power-law behaviour, with some evidence for a mass 
dependence of the clustering amplitude, in particular at the up- 
per mass limit. The overall shape is well described by a rather 
flat power-law ^(s) ~ (so/sy, with slope y ~ I and a correla- 
tion length sq between 6 and ~ 10 h ' Mpc. These values for 
the shape and amplitude of ^(s) are similar to those measured 
for luminous red galaxies at z = 0.55 in the 2SLAQ survey (see 
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Fig. 20. The redshift-space, angle-averaged correlation function 
^(s) of the mass-selected samples in the redshift bin [0.5, 1]. A 
mild systematic mass dependence is visible. 



Fig. 7 in lRoss et al .l (2001)) and for lum inous early-type galax- 
ies m the 2dFGRS (Norberg etalllZQOj. This is consistent with 
the most massive objects in the lOK sample being predominantly 
red, early-type galaxies which show moderate or no evolution in 
the overall clustering amplitude with redshift. 

In Fig.|2Tl instead, we compare ^{s) of our "reference" sam- 
ple with log(M/MQ) > 10, with the mean and scatter (at Icr 
and 3cr confidence, respectively) of the similarly-selected set 
of BDL40 mocks. Despite the angular average, we note a be- 
haviour which is similar to that observed in Wp{rp) although 
now the agreement extends to slightly larger scales. The ob- 
served clustering is compatible with the predictions of the stan- 
dard model (to better than the 68% level) on scales smaller than 
~ 2 h ' Mpc. On larger scales, also ^(s) shows excess power 
with respect to the models, which places the zCOSMOS volume 
at the upper 3(t limit of the statistical distribution obtained from 
the mocks. This exercise shows that even after angle-averaging 
our clustering estimator, the amount of structure present in this 
specific volume of the Universe remains outstanding in com- 
parison to the model expectations. The conclusion can only be 
that either we have been very unlucky in the selection of the 
COSMOS field and picked up a fluctuation which has a prob- 
ability of ~ 1% to be found in such a volume, or fluctuations 
with this amplitude are in reality more common than what the 
standard cosmology predicts. 

12. Summary 

We have used the lOK zCOSMOS spectroscopic sample to study 
galaxy clustering as a function of galaxy luminosity and stellar 
mass, in the range of redshift [0.2,1]. 

To this end, we built luminosity and mass-selected sam- 
ples from the lOK catalogue sampling three separate redshift 
ranges. We used mock catalogues to quantify the effect of stellar 
mass incompleteness on the measured clustering, as a function 
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s [Mpc/h] 

Fig. 21. Comparison of the redshift-space correlation function 
for galaxies with log(M/M0) > 10 at 0.5 < z < 1 in the zCOS- 
MOS sample (filled circles), with the model predictions from the 
BDL40 mock samples. The solid line gives the average of the 40 
mocks, with the dashed areas corresponding to the Icr and 3cr 
confidence error corridor Similarly to what we found for Wp(rp), 
the agreement between the models and the zCOSMOS measure- 
ment at 1 h"' Mpc is remarkable. Nevertheless, the large-scale 
shape of the zCOSMOS ^(s) is at ~ 3cr the mean amplitude of 
the mock catalogues. 

of redshift. We carefully checked our covariance and error esti- 
mate techniques, comparing the performances of methods based 
on the scatter in the mocks and on bootstrapping schemes. We 
adopted the latter, based on 200 resamplings of 8 sub-volumes 
of the survey, as the most appropriate description of of the co- 
variance properties of the data. 

By measuring the redshift-space correlation functions ^(s) 
and ^(rp,7r) and the projected function Wp{rp) for these sub- 
samples, we found the following results. 

- Surprisingly, we do not see any clear dependence on lumi- 
nosity of the correlation function at all redshifts. This is 
at odd with results in the local Universe by the 2dFGRS 
and with mesurements at similar redshift by the VVDS and 
DEEP2 surveys, that found a significant steepening of Wp(rp) 
with luminosity. 

- We find a mildly more evident (although not striking) depen- 
dence of Wp(rp) on stellar mass, especially on small scales. 
The central redshift bin (0.5 < z < 0.8), displays in general a 
more evident effect, with a very flat shape of Wp{rp) on scales 
rp - [1 - 10] h^^ Mpc. The overall shape of the correspond- 
ing map of ^{rp, n) shows strong distortions that we interpret 
as the effect of dominant structure extending preferentially 
perpendicularly to the line of sight. 

- From comparison to the SDSS measurements at z ~ 0, 
we do not see any significant evolution with redshift of the 
amplitude of clustering for bright and/or massive galaxies. 
Together with previous results from VVDS, this is consis- 
tent with a more rapid evolution of the linear bias for the 



most massive objects, with respect to the general population. 
In the zCOSMOS sample this invariance in the clustering 
amplitude between z ~ 1 and z ~ seems to remain valid 
down to smaller masses than in the VVDS, an effect easily 
explained by the overall larger clustering amplitude observed 
in general for z > 0.5 in this sample. This is evidenced by 
a much flatter shape (higher amplitude) of Wpifp) of zCOS- 
MOS galaxies with respect to VVDS galaxies, when selected 
with the same criteria. 
- This particularly high level of structure is confirmed by com- 
parison of the measured Wp{rp) and ^{s) at 0.5 < z < 1 
with model predictions, concentrating on the sample with 
log(M/M0) > 10. On scales smaller than ~ 1 -2 h"' Mpc, the 
observations agree very well with the model expectation val- 
ues in the standard ACDM scenario for a linar bias b ~ 1.6. 
On these scales, the measured values are compatible to bet- 
ter than 68% with the BDL40 mocks. On larger scales, how- 
ever, the observed clustering amplitude is reproduced in only 
a few percent of the mocks. In other words, if the shape of the 
power spectrum is that of ACDM and the bias has no "innat- 
ural" scale-dependence, COSMOS has picked up a volume 
of the Universe which is rare, 2 - 3a- positive fluctuation. 
This conclusion is corroborated also by comparison with the 
VVDS measurements, which on the other hand lie on the 
lower side of the distribution, at about 1.5- Icr. 

Acknowledgements. We thank the anonymous referee for a detailed review 
of the manuscript that helped to improve the paper We thank G. De Lucia, 
J. Blaizot, S. Phleps and A. Sanchez for their thorough comments on an early ver- 
sion of the manuscript. This work was supported by grant AS1/COFIS/WP3110 
1/026/07/0. 



References 

Abbas, U., & Sheth, R. K. 2005, MNRAS, 364, 1327 

Abbas, U., & Sheth, R. K. 2006, MNRAS, 372, 1749 

Abbas, U., & Sheth, R. K. 2007, MNRAS, 378, 641 

Abbas, U., de la Torre, S., Le Fevre, O. et al. 2009, A&A, submitted 

Abbas, U. et al. 2009, A&A, in preparation 

Arnouts, S., Walcher, C. J., Le Fevre, O. et al. 2007, A&A, 476, 137 

Bell, E. F, Naab, T., Mcintosh, D. H. et al. 2006, ApJ, 640, 241 

Benoist, C, Maurogordato, S., da Costa, L. N., Cappi, A., & Schaeffer, R. 1996, 

ApJ, 472, 452 
Bernstein, G. M., 1994, ApJ, 424, 569 

Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 476 
Bolzonella, M. et al. 2009, A&A, in preparation 

Blaizot, J., Wadadekar, Y, Guiderdoni, B. et al. 2005, MNRAS, 360, 159 
Brown, M. J. 1., Zheng, Z., White, M. et al. 2008, ApJ, 682, 937 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Capak, R, Aussel, H., Ajiki, M. et al. 2007, ApJS, 172, 99 
Chabrier, G. 2003, PASP, 115, 763 

Coil, A. L., Newman, J. A., Cooper, M. C, Davis, M., Faber, S. M., Koo, D. C, 

& Willmer, C. N. A. 2006, ApJ, 644, 671 
Coil, A. L., Newman, J. A., Croton, D. et al. 2008, ApJ, 672, 153 
Cooray, A., & Sheth, R. 2002, Phys. Rep., 372, 1 
Conroy, C, Wechsler, R. H., & Ki'avtsov, A. V. 2006, ApJ, 647, 201 
Daddi, E., Rottgering, H. J. A., Labbe, 1. et al. 2003, ApJ, 588, 50 
Davis, M., & Geller, M. J. 1976, ApJ, 208, 13 
Davis, M., & Peebles, J. P E. 1983, ApJ, 267, 465 
de la Torre, S., Le Fevre, O. et al. 2009, in preparation 
De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 

Feldmann, R., CaroUo, C. M., Porciani, C. et al. 2006, MNRAS, 372, 565 
Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. 1994, MNRAS, 
266, 50 

Giacconi, R., Zirm, A., Wang, J. et al. 2002, ApJS, 139, 369 
Giovanelli, R., Haynes, M. P, & Chincarini, G. L. 1986, ApJ, 300, 77 
GarilU, B., Le Fevre, O., Guzzo, L. et al. 2008, A&A, 486, 683 
Guzzo, L., Strauss, M. A., Fisher, K. B., GiovanelU, R., & Haynes, M. P 1997, 
ApJ, 489, 37 

Guzzo, L., Bardett, J. G., Cappi, A. et al. 2000, A&A, 355, 1 
Guzzo, L., Cassata, P, Finoguenov, A. et al. 2007, ApJS, 172, 254 
Guzzo, L., Pierleoni, M., Meneux, B. et al. 2008, Nature, 451, 541 



18 



Meneux, B., et al.: Clustering as a function of galaxy luminosity and mass 



Hamilton, A. J. S. 1988, ApJ, 331, L59 

lovino. A., Giovanelli, R., Haynes, M., Chincaiini, G., & Guzzo, L. 1993, 

MNRAS,265. 21 
Ilbert, O., Tresse, L., Zucca, E. et al. 2005, A&A, 439, 863 
Ilbert, O., Capak, R, Salvato, M. et al. 2009, ApJ, 690, 1236 
Kaiser, N., 1987, MNRAS, 227, 1 

Kitzbichler, M. G., & White, S. D. M. 2007, MNRAS, 376, 2 
Kovac, K., Lilly, S. J., Cucciati, O. et al. 2009, arXiv: 0903.3409l 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 

Le Fevre, O., Saisse, M., Mancini, D. et al. 2003, Proc. SPIE, 4841, 1670 

Le Fevre, O., Vettolani, G., Gaiilli, B. et al. 2005a, A&A, 439, 845 

Le Fevre, O., Guzzo, L., Meneux, B. et al. 2005b, A&A, 439, 877 

Li, C, Kauffmann, G., Jing, Y. R, White, S. D. M.. Borner G., & Cheng, F. Z. 

2006, MNRAS, 368,21 
Lilly, S. J., Le Fevre, O., Renzini, A. et al. 2007, ApJS, 172, 70 
Lilly, S. J., Le Brun, V., Maier, C. et al. 2009, ApJS, submitted 
Ling, E. N., Ban-ow, J. D., & Frenk, C. S. 1986, MNRAS, 223, 21L 
Maurogordato, S., & Lachieze-Rey, M. 1991, ApJ, 369, 30 
McCracken, H. J., Peacock, J. A., Guzzo, L. et al. 2007, ApJS, 172, 314 
McCracken, H. J., Ilbert, O., Mellier, Y. et al. 2008, A&A, 479, 321 
Meneux, B., Le Fevre, O., Guzzo, L. et al. 2006, A&A, 452, 387 
Meneux, B., Guzzo, L., Garilh, B. et al. 2008, A&A, 478, 299 
Norberg, R, Baugh, C. M., Hawkins, E. et al. 2001, MNRAS, 328, 64 
Norberg, R, Baugh, C. M., Hawkins, E. et al. 2002, MNRAS, 332, 827 
Norberg, R, Baugh, C. M., Gaztanaga, E., & Croton, D. J. 2008, 

arXiv:0810.1885 

Peebles, P. J. E. 1980, The Large Scale Structure of the Universe (Princeton 

University Press, Princeton) 
Phleps, S., Peacock, J. A., Meisenheimer, K., & Wolf, C, 2006, A&A, 457, 145 
Polio, A., Meneux, B., Guzzo, L. et al. 2005, A&A, 439, 887 
Polio, A., Guzzo, L., Le Fevre, O. et al. 2006, A&A, 451, 409 
Porciani, C, & Giavalisco, M. 2002, ApJ, 565, 24 
Porciani, C. et al. 2009, ApJ, in preparation 

Pozzetti, L., Bolzonella, M., Lamareille, F. et al. 2007, A&A, 474, 443 
Pozzetti. L. et al. 2009. A&A, in preparation 
Press, W. H., & Schechter R 1974, ApJ, 187, 425 
Rettura, A., Rosati, R, Strazzullo, V. et al. 2006, A&A, 458, 717 
Ross, N. R, da Angela, J., Shanks, T. et al. 2007. MNRAS, 381, 573 
Scoccimarro, R. 2004, Phys. Rev. D, 70, 083007 
Scoville, N., Abraham, R. G., Aussel, H. et al. 2007, ApJS, 172, 38 
Scoville, N., Aussel, H., Benson, A. et al. 2007, ApJS, 172, 150 
Skibba, R. A., Sheth, R. K., & Martino, M. C. 2007, MNRAS, 382, 1940 
Smith, R. E., Peacock, J. A., Jenkins, A. et al. 2003, MNRAS, 341, 1311 
Springel, V., White, S. D. M., Jenkins, A. et al. 2005, Nature, 435, 629 
Springel, V., Frenk, C. S., & White, S. D. M. 2006, Nature, 440, 1137 
Stringer M. J., Benson, A. J., Bundy, K., Ellis, R. S., & Quetin, E. L. 2009, 
MNRAS, 143 

Tinker, J. L., Norberg, R, Weinberg, D. H., & Wairen, M. S. 2007, ApJ, 659, 877 

Wake, D. A., Sheth, R. K.; Nichol, R. C. et al. 2008, MNRAS, 387, 1045 

Wang, L., Li, C, Kaulfmann, G.. & De Lucia, G. 2006, MNRAS, 371. 537 

Wang, L., Li, C, Kauffmann, G., & De Lucia, G. 2007, MNRAS, 377, 1419 

Wan-en, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 2006, ApJ, 646, 881 

Willmer, C. N. A., da Costa, L. N., & Pellegrini, R S. 1998, AJ, 115, 869 

White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 

Zehavi, I., Blanton, M.R., Frieman, J. A. et al. 2002, ApJ, 571, 172 

Zehavi, I., Zheng, Z., Weinberg, D. H. et al. 2005, ApJ, 630, 1 

Zheng, Z., Coil, A. & Zehavi, I. 2007, ApJ, 667, 760 

Zucca, E., Ilbert, O., Bardelli, S. et al. 2006, A&A, 455, 879 

Zucca, E., Bardelli, S.. Bolzonella, M. et al. 2009. A&A, in preparation 

' Max-Planck-lnstitut fiir Extraterrestrische Physik, 
Giessenbachstrasse, 85748 Garching-bei-Miinchen, Germany 

^ Universitats-Stemwarte, Scheinerstrasse 1, Munich D-81679, 
Germany 

' INAF-Osservatorio Astronomico di Brera, Via Bianchi 46, I- 
23807 Merate (LC), Italy 

Laboratoire d'Astrophysique de Marseille, UMR 6110 CNRS 
Universite de Provence, BPS, F-13376 Marseille Cedex 12, France 

' INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica, Via 
Bassini 15, 1-20133 Milano, Italy 

* Argelander Institute for Astronomy, Auf dem Hiigel 71, D-53121, 
Germany 

' Institute of Astronomy, ETH Zurich, Zurich, Switzerland 
^ INAF - Osservatorio Astronomico di Bologna, Bologna, Italy 
' INAF - Osservatorio Astronomico di Torino, Strada Osservatorio 
20, 1-10025 Pino Torinese (TO), Italy 



INAF - Osservatorio Astronomico di Brera, via Brera 28, Milano, 
Italy 

" Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de 

Toulouse, CNRS Toulouse, F-31400, France 

European Southern Observatory, Garching, Germany 
Dipartimento di Astronomia, Universita di Padova, Padova, Italy 
Dipartimento di Astronomia, Universita degli Studi di Bologna, 

Bologna, Italy 

Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA 

Berkeley Lab & Berkeley Center for Cosmological Physics, 
University of California, Berkeley, CA 94720, USA 
" Centre de Physique Theorique, UMR 6207 CNRS Universite de 
Provence, F- 13288 Marseille, France 

Institut d'Astrophysique de Paris, Universite Pierre & Marie 
Curie, Paris, France 

" INAF - Osservatorio Astronomico di Roma, via di Frascatti 33, 
00040 Monte Porzio Catone, Italy 



